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Abstract: In nuclear reactors, temperature fluctuations of fluids may cause fatigue 
damage to adjacent structures; this is referred to as thermal striping. Research on 
thermal striping in the upper plenum has mainly focused on fluid fields. Few 
experimental studies have been reported on solid structures in a fluid field with a 
coaxial jet. This study entailed an experimental study of the temperature fluctuations in 
the fluid and on a plate surface caused by a coaxial jet. The temperature fluctuations of 
the fluid and plate surfaces located at different heights were analyzed. The cause of the 
temperature fluctuation was analyzed using a transient temperature distribution. The 
results show that the mixing of the hot and cold fluids gradually becomes uniform in 
the positive axial direction. The average surface temperatures tended to be consistent. 

When the jet reaches the plate surface, the swing of the jet center, contraction and 
expansion of the cold jet, and changes in the jet shape result in temperature fluctuations. 
The intensity of the temperature fluctuation was affected by the position. More attention 
should be paid when the plate is located at a lower height, and between the hot and 
cold-fluid nozzles. 


Keywords: temperature fluctuation, thermal striping, coaxial jet, thermal mixing, the 
upper plenum of nuclear reactor 


1. Introduction 


During the operation of nuclear reactors, owing to the uneven temperature of fluids, 
mixing of fluid streams with different temperatures, alternating thermal stress, and 
severe fatigue damage may occur in adjacent structures |'7!, This phenomenon, referred 
to as thermal striping, is pronounced in the upper plenum of nuclear reactors'*!, pipe 
elbows |“, and T-junction P! Elbow failure is often caused by fluid stratification at 
different temperatures'®!, The causes of thermal striping at the T-junction and upper 
plenum of nuclear reactors are similar, namely, the mixing of fluids at different 
temperatures!7*!, Because temperature change is the main cause of damage, the analysis 
of temperature fluctuations at the aforementioned locations is the basis for the study of 
thermal striping ”. 

Because T-junctions are a common structure in pipelines, temperature fluctuations 
in these locations have been extensively studied. At T-junctions, two fluids are mixed 
vertically at different temperatures. Temperature fluctuations occur downstream of the 


mixing position!!°!, Several studies have been reported on temperature fluctuations in 


T-junctions. Most such studies have focused on the fluid temperature distribution inside 
T-junctions under different conditions. Kimura et al. !''! characterized the jet 
impingement and mixing between a hot stream and a smaller cold branch flow. It was 
found that the temperature fluctuation in the T-junction with the elbow upstream had a 
large component at a low frequency compared to the straight case. Zhou et al.!!?] 
conducted thermal mixing experiments to investigate the thermal striping close to a 
weld connection downstream of a 90° T-junction. The results show that the change in 
temperature difference has a significant influence on the stability of the thermal 
stratification in the mixing flow, as well as on the phenomenon of reverse flow. Walker 
et al.!!3! and Kickhofel et al.!'4! characterized T-junction mixing by using wire mesh 
sensors. The experiments obtained a large amount of temperature distribution data on 
the cross-section of the T-junction; these data are ideal for verifying the accuracy of the 
simulation. On the basis of experimental results, many scholars have performed 
simulation analyses of temperature fluctuations in T-junctions. Most studies indicate 
that large eddy simulations (LES) can reproduce long-term temperature fluctuations at 
T-junctions [!5' 16], Qian and Kasahara! "l performed LES analyses of flow and 
temperature fields at a T-junction by using different subgrid-scale (SGS) models. The 
simulation results indicate that the numerical schemes significantly affect the 
temperature distribution and temperature-fluctuation intensity. An approach using the 
dynamic Smagorinsky model (DSM) and a hybrid scheme with a large blending factor 
can accurately predict fluid temperature fluctuations. In addition, some studies have 
focused on the heat transfer between fluids and structures. Guo et al.!'®! measured the 
outer wall temperature and solved one-dimensional and two-dimensional inverse heat 
conduction problems (IHCPs) to estimate the temperature fluctuations on the inner wall. 
LES can also be used for the conjugate heat transfer (CHT) analysis between the fluid 
and pipe wall in a T-junction |!?!, Utanohara et al.!?°! performed a fluid—structure 
coupled simulation of a CHT to investigate the predictive performance of the flow and 
temperature fields and the temperature fluctuation on the inner surface of the pipe at 
the T-junction. The simulation reproduced the trend of the experimental data, 
particularly the velocity fluctuation intensity peak near the wall. Although the mixing 
form of the T-junction differs from that of the upper plenum, these studies provide 
important methodological references for studying the temperature fluctuations in the 
upper plenum. 

In the upper plenum, owing to the distribution characteristics of the core 
assemblies, there were many fuel assemblies around each control-rod assembly at 
different distances. Moreover, the coolant flow in each channel around the fuel 
assembly is not always consistent. This led to different coolant temperature profiles in 
each channel, resulting in thermal striping. In the upper plenum, multiple fluids of 
different temperatures flow and mix in parallel, causing temperature fluctuations in the 
flow direction! ?! l. Thermal striping in the upper plenum damages many vital 
components. Therefore, they have attracted significant attention!*”! Over the past few 
decades, various thermal striping studies focusing on the upper plenum of nuclear 
reactors have been reported. The upper plenum was previously simplified as a parallel- 


jet model. Wakamatsu et al.'?3! conducted an experimental study on thermal striping in 
an upper plenum. In their study, a parallel-jet model was adopted to analyze the 
attenuation characteristics of temperature fluctuations with water and liquid sodium, 
and a general attenuation range was provided. Tokuhiro et al.!'! further investigated 
thermal striping by using a triple-parallel jet model. The mixing locations were 
measured and determined using the same velocities and various velocity cases. Chacko 
et al.P® conducted a numerical study by using a triple parallel jet. It was found that LES 
could be used to analyze the unsteady characteristics of thermal striping. Their study 
also indicated the limitations of the Reynolds-averaged Navier-Stokes (RANS) 
approach for unsteady heat transfer simulations. Wang et al.!?>°! used LES to study the 
temperature fluctuations caused by a parallel jet in different media. They concluded that 
the temperature-fluctuation amplitudes of the lead-based materials were more 
significant than those of sodium. Yu et al. [?7] conducted steady and unsteady 
simulations (RANS, unsteady RANS (URANS), and LES) based on the triple-parallel 
jet model. The temperature fluctuation in the unsteady simulation agreed well with the 
experimental results. Lomperski et al.!78! presented velocity and temperature field 
measurements in which two air jets were mixed and impinged on a wall. This indicates 
that the mixing is better when the jet impacts the wall than when it is freely extended. 
Krishna et al.'??! employed two parallel jets impinging on a lattice plate above them. 
The Reynolds stress turbulence model was used to evaluate temperature fluctuations 
near the plate. Jets with a unity velocity ratio exhibited the maximum temperature 
fluctuations. Kimura et al.°°! employed a triple parallel jet to experimentally investigate 
the transfer characteristics of temperature fluctuation from the fluid to the structure. 
The nonstationary heat transfer characteristics can be represented by a heat transfer 
coefficient that is constant in time and independent of the frequency of the temperature 
fluctuation. The aforedescribed studies analyzed the thermal striping phenomenon by 
using parallel-jet and triple-parallel jet models, which can reflect the two-dimensional 
characteristics of the temperature fluctuations. However, the thermal mixing in the 
upper plenum is a complex phenomenon. Thermal striping exhibits strong three- 
dimensional characteristics at this location. 

Consequently, many studies have simplified the upper plenum of nuclear reactors 
to a coaxial-jet model. Because coaxial jets can reflect the three-dimensional 
characteristics of temperature fluctuations, their use to describe the thermal striping 
phenomenon is reasonable!?!!, Lu et al.!**! conducted coaxial-jet experiments to 
investigate the temperature-fluctuation characteristics of a flow field. They found that 
the alternating extrusion of hot- and cold-flow columns caused temperature fluctuations. 
Kok et al.'?! experimentally and numerically investigated the mixing behavior of hot 
and cold fluids for coaxial-jet flows. The thermal mixing efficiency increased with the 
temperature difference between the hot and cold jet flows. The thermal mixing 
performance was improved by increasing the hot-jet flow rate. Tenchine et al.!*+! 
conducted jet experiments on sodium and air. The objective of their experiment was to 
evaluate the possibilities and limitations of simulating sodium-mixing jets by using a 
typical fluid. It was found that a rather good estimation of the mean sodium temperature 
and temperature fluctuation was possible when tests using a simulant fluid were 


performed at a Reynolds number close to that of the prototype. Their study on the 
temperature fluctuation induced by a coaxial jet was limited to the fluid field. The swing 
of the jet likely causes a constant change in the fluid temperature at some positions. 
However, this change is also transferred to the adjacent solid structures, resulting in 
temperature fluctuations in solid structures. However, the plate in the fluid causes the 
flow field to change. Therefore, to elucidate the mechanism of thermal striping, the 
temperature fluctuations of the flow field and structure surface caused by coaxial-jet 
flows with an adjacent solid structure must be experimentally examined. 

In this study, the thermal striping phenomenon induced by a coaxial jet was 
experimentally investigated using a coaxial-jet experimental setup. The temperature 
fluctuations of the fluid and plate surfaces in the coaxial-jet flow field were examined. 
The transient temperatures of the plates located at different heights were measured and 
analyzed. Subsequently, variations in the average temperature and local temperature- 
fluctuation intensity along the axial direction were studied. 


2. Experimental setup 


2.1 System description 

An experimental setup based on a coaxial-jet model was designed and built to 
study the temperature fluctuations. The setup included hot fluid, cold fluid, mixing, and 
measurement systems, as shown in Fig. l(a). Both hot- and cold-fluid systems consist 
of a thermostatic water tank, booster pump, electromagnetic flowmeter, storage tank, 
valves, and pipes, which can provide steady water flow at specified temperatures and 
flow rates. The mixing system included coaxial-jet nozzles, a mixing zone, valves, and 
pipes. The coaxial-jet nozzle was designed using two independent passages. Cold fluid 
flows through the central passage, whereas hot fluid flows through the outer passage. 
The two fluids are mixed coaxially in the mixing zone. The measurement system is 
mainly used for the temperature measurement of the fluid and plate surface in the 
mixing zone and comprises thermocouples, brackets, a data acquisition card, and a 
computer. 

During the experiment, hot and cold water were pumped into the mixing zone from 
the thermostatic water tank via the valves. The mixed water flows into the storage tank 
and then recycled into the thermostatic water tank. Fig. 1(b) shows a photograph of the 
experimental setup. 


2.2 Test zone and inlet temperature test point 

The test zone was a transparent acrylic box with dimensions 250 mm x 250 mm x 
100 mm, as shown in Fig. 1 (c). The coaxial jet inlet is located at the bottom of the test 
zone. The nozzle was placed 10 mm above the bottom to reduce the impact of the 
bottom plate on the inlet fluid. The outlet was located on the two sides above the box, 
away from the test zone, to reduce its impact on the flow field. The brackets were 
movable and can be placed at certain heights within the box to enable temperature 
measurements at different heights. 

The cold-fluid inlet was circular with a diameter of 10 mm. The hot-fluid inlets 


were annular with 24.5 mm and 32 mm. Three measurement points were evenly 
arranged on the circumferences of the hot and cold-fluid nozzles to measure the inlet 
temperatures of the fluids, as shown in Fig. 1(d). 
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Fig. 1 Experimental setup and measurement points 


2.3 Uncertainty analysis 

The temperature was measured using a T-type thermocouple with a diameter of 
0.08 mm. The time interval for temperature sampling was 0.01 s. The velocity was 
calculated using the hot and cold flow rates and the cross-sectional area of the coaxial- 
jet nozzle. The hot and cold flow rates were measured using electromagnetic flow 
meters. The setup temperature measurement range was 20-90 °C, and the flow rate was 
0.141—2.802 m°/h. The temperature measurement uncertainty was determined using a 
thermocouple and a data acquisition card. An electromagnetic flow meter was used to 
directly determine the uncertainty in the flow measurement. 

The measurement range and accuracy of the equipment were as follows. The 
measurement range of the T-type thermocouple was 0-150 °C, and the expanded 
uncertainty was U(Tihermocouple) = 0.2 K with an inclusion factor k = 2. The standard 
uncertainty of the thermocouple u(Tinermocouple) = 0.1 K. The interval half-width of the 
data acquisition card (model, Hioki) a(Tyioxi) = 0.5 K for the T-type thermocouple 
within 0-100 °C. The electromagnetic flowmeter measurement range was 0-3 m°/h 


with a half-width a(Viowmerer) = 0.0102 m3/h. Considering uniform distribution, the 
inclusion factor was taken as k=/3, and the standard uncertainty of the data acquisition 
card u(THioki) = 0.28 K. The standard uncertainty of flowmeter u(Vfowmeter) = 0.00589 
m°/h. The combined standard uncertainty of the temperature measurement u,(7) can 
be obtained using Eq. (1) B5, 


u-(T) 3 | E A a +u(T hior) (1) 


Thus, the combined standard uncertainty of the temperature u.(7) =0.311 K. Similarly, 
the standard uncertainty u(V) and relative expanded uncertainty range u,,.;,(V) of the 
flow measurement can be obtained using Eqs. (2) and (3), respectively: 


u(V) = u(Vjtowmeter) (2) 
u(V) uV) 
Ucrels(V) = k( Pox a Vpis ) (3) 
The standard uncertainty u(V) = 0.00589 m?/h, and the expanded uncertainties of the 
flow measurement Uecrerl V) = 0.42%-8.4%. 
The inlet temperature T; is the average of the values at the three measurement 
points. 
The combined standard uncertainty of the inlet temperature measurement u,(T7;,) can 
be calculated using Eq. (4). 


2 2 
u(Tin) = EE pero TT iT (4) 


The combined uncertainty of inlet temperature measurement u,(7;,,) was 0.180 K. 
This shows that the setup has high measurement stability and can be used for 
experimental measurements. 


2.4 Data processing method 
For the convenience of discussion, the temperature and radius values were 


nondimensionalized. The dimensionless temperature is expressed as 
«FT coi 
at ee sae (5) 
T Hot™ T Cold 
where Tcoa and Typos denote the inlet temperatures of the cold and hot fluids, 
respectively. The average values measured at the test points at the inlet were used. The 


dimensionless average temperature was calculated as 


n 
* 1 * 
Tiyg== )T 6) 
i=l 
The dimensionless radius is the ratio of the diameter to the cold-fluid inlet diameter. 
Pe 2 
r= 7 
’Cold a) 


where /cojq is the radius of the cold-fluid inlet (5 mm). The root-mean-square (RMS) 
temperature can be used to express the intensity of temperature fluctuations. 


(8) 


3. Coaxial-jet experiment scheme 


The temperature sensors were fixed in the fluid and on the plate surface using 
detachable brackets. The temperature measurements of the plate surface were divided 
into two parts. First, the 29 points shown in Fig. 1(e) (A) were arranged to measure the 
entire plane of the lower surface. Because the temperature fluctuation caused by the 
coaxial jet was symmetrical, the measurement points were arranged along the radial 
direction in further tests. There were 13 measurement positions on the plate surface, as 
shown in Fig. 1(e)(B), at which the thermocouples were pasted without covering so that 
the thermocouple can be in contact with the fluid. To measure the fluid temperature, 13 
thermocouples were fixed on a tungsten steel wire under the plate surface, as shown in 
Fig. 1(e) (C). The measurement point positions were the same as those shown in Fig. 
1(e) (B). 

The plate size was 100 mmx100 mm, and the thickness was 1 mm. Its material is 
stainless steel 316, with a density of 7.98 g/cm?, thermal conductivity of 15.1 W/(m:K), 
and specific heat capacity of 0.502 J/(g:K). 

For all the test cases, the flow velocity was 0.5 m/s. The plate height z was varied 
from 30 to 110 mm. The cold and hot water temperatures were set at 20 °C and 40 °C, 
respectively. The thermocouple sampling time interval was 0.01 s, and the measurement 
time of a single test was approximately 180 s. The experimental conditions are listed in 
Table 1. 


Table 1 Experiment conditions 


Heich Hot water Cold water 
aay Temperature Velocity Temperature Velocity 
(°C) (m/s) (°C) (m/s) 
Casel 30 40 0.5 20 0.5 
Plate Case2 40 40 0.5 20 0.5 
surface Case3 50 40 0.5 20 0.5 
Case4 60 40 0.5 20 0.5 
. Case5 10 40 0.5 20 0.5 
Fme ae 2 40 0.5 20 0.5 
in radial 
direction Case7 30 40 0.5 20 0.5 
Caseg 40 40 0.5 20 0.5 
Plate Case9 30 40 0.5 20 0.5 
pene Casel0 40 40 0.5 20 0.5 
in the Casel1 50 40 0.5 20 0.5 
dial Case12 70 40 0.5 20 0.5 
n adiar Casel3 90 40 0.5 20 0.5 
ection Casel4 -110 40 0.5 20 0.5 


4. Temperature fluctuation on the plate surface 


Considering a plate height of 50 mm as an example, in addition to the temperature 
fluctuations on the plate surface, the fluid temperature fluctuation between the coaxial- 


jet nozzle and the plate was analyzed. The height measurement range for the fluid 
temperature was 10—40 mm, shown as by Z1—Z4 in Fig. 1(e). The dimensionless heights 
were calculated as 
« 2f 
Jas ©) 
where zy is the measured height of the fluid. zp is the plate height, which was 50 mm in 


this study. 


4.1 Transient temperature analysis 

Subsequent average temperature analysis revealed that the dimensionless average 
temperature of the fluid was symmetric. Therefore, an analysis of the transient 
temperature was performed along the radius (P15—P18 in Fig. 1(e) (A)). The 
temperature fluctuations were entirely random, as shown in Fig. 2(a). When the plate 
height zp = 50 mm, the temperature at each radial position significantly fluctuated; r°= 
0 showed the most pronounced fluctuations. In addition, there were small fluctuations 
overall temperature increase and decrease. For example, at r* = 0, the temperature 
shows a decreasing trend between 0 and 0.2 s; however, this process is accompanied by 
two minor temperature increases. 

The cause of the temperature fluctuation can be analyzed by comparing the 
temperature fluctuation in Fig. 2(a) with the transient temperature distribution on the 
plate surface at different times shown in Fig. 3(c). These temperature fluctuations can 
mainly be attributed to the fluid instability. As indicated by the transient temperature 
distribution on the plate structure, this instability is mainly reflected in the following 
three aspects. First, the cold fluid would expand and contract as time progresses. For 
example, the cold fluid expands during 0-0.2 s and contracts during 0.2-0.4 s. 
Correspondingly, the center (r* = 0) temperature decreases in 0-0.2 s, and the 
temperature increases in 0.2—0.4 s (Fig. 2(a)). Second, the jet center would swing over 
time. At 0.2 s, the center of the cold fluid was at r“ = 0, as shown in Fig. 3(c). At 0.5 s, 
the center moved downward. At 0.6-0.7 s, the cold-fluid center appeared on the left. 
The transient temperatures are presented in Fig. 2(a). The temperature at r= 2 increases 
when the cold-fluid center swings to the left at 0.4 s. Finally, the shape of the jet also 
changes significantly. Occasionally, the temperature distribution on the plate surface 
presents a relatively regular circle, for instance, that at 1.1 s. It may also extend in 
different directions, for example, the observations corresponding to 0.1 s, 0.2 s, and 0.9 
s. Mostly, the shape of the jet is relatively irregular. Under the combined influence of 
these three factors, the temperature would fluctuate at each point on the plate surface. 

In the flow field between the jet nozzle and the plate, the transient temperatures at 
zy = 0.2 and zr = 0.6 were as shown in Fig. 2(b). The position at which the temperature 
fluctuation is relatively prominent changes along the axial direction. The position with 
notably evident temperature fluctuation is at r*= 1 when zy = 0.2 and at r’ = 0 when zy 
= 0.6. Compared to plate surface, temperature changes were more frequent in the flow 
field zy = 0.2 and zř = 0.6. 
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Fig. 2 Dimensionless transient temperature, average temperature, and intensity of temperature 
fluctuations of fluid on the plate surface 


4.2 Average temperature distribution 

Over a long period, the temperatures at various radial positions fluctuated around 
a fixed value. This is reflected in the dimensionless average temperature values. The 
dimensionless average temperature of the plate surface at zp = 50 mm is shown in Fig. 
5(a). The average temperature was the lowest at the center and continued to increase 
along the radial direction. The temperature distribution was almost circular. Therefore, 
the temperature along the radial direction was selected for further analysis. The 


dimensionless average temperature of the fluid in the radial direction was symmetrical, 
as shown in Fig. 2(c). The blue and red blocks on the horizontal axis represent the 
positions of the cold- and hot-fluid nozzles, respectively. The average temperature was 
the lowest at the center. In the positive radial direction, it increased rapidly and reached 
the maximum value above the hot-fluid nozzle. Subsequently, the average temperature 
decreased slowly and finally reached a constant value of approximately 0.836. 
Moreover, the lower the height, the smaller would be the minimum value and the larger 
would be the maximum value. The average dimensionless temperature on the plate 
surface showed a considerably different behavior. It increased continuously In the 
positive radial direction. There was no longer a prominent high-temperature peak above 
the hot-fluid nozzle. This phenomenon may be attributed to the cold fluid in the center 
flowing in all directions and rapidly mixing with the hot fluid when the jet reaches the 
plate surface. 

To examine the temperature variation in the axial direction, the center position 
(P15 in Fig. 1(e)(B), r“ = 0), mixing zone between the hot and cold fluids (P16 in Fig. 
1(e)(B), r= 2), and outer position (P18 in Fig. 1(e)(B), r*= 6) were considered for the 
analysis. Fig. 2(d) shows the average dimensionless temperature of the fluid in the axial 
direction. The average temperature increased at the center (r= 0) in the axial direction. 
It decreases at r= 2 and remains constant at r= 6. This variation is linear between zy 
= 0.2-0.8. 


4.3 Intensity of temperature fluctuation 

Intensity of the temperature fluctuations can better characterize their severity. This 
can be expressed as the RMS temperature. For zp = 50 mm, the RMS dimensionless 
temperature on the plate surface was as shown in Fig. 5(b). The temperature evidently 
fluctuated above the jet nozzle. The RMS temperature along the radial direction is 
shown in Fig. 2(e). At zr = 0.2, 0.4, and 0.8, the RMS temperature initially increased 
and then decreased, with the maximum value observed at r*= 1. At z = 0.6, the RMS 
temperature decreased continuously from r* = 0. At the plate surface, the RMS 
temperature was relatively higher in the range of r*< 4. 

Fig. 2(f) shows the variation of the intensity of the temperature fluctuations in the 
axial direction. At r’ = 0, the RMS temperature increased when zf < 0.6 and gradually 
decreases when zf > 0.6. It shows an increasing trend at r*= 2 and remains constant at 
r = 6. Severe temperature fluctuations indicated that the hot and cold fluids constantly 
alternated at this position. This point was the junction of the hot and cold fluids. 
Therefore, the intensity of the temperature fluctuation at r*= 0 reflected the amplitude 
of the jet swing. When zy < 0.6, the jet swing amplitude increased in the axial direction. 
The junction between the hot and cold fluids can swing toward the center. Therefore, 
the RMS temperature at r*= 0 continued to increase. When zy > 0.6, the plate makes 
the jet more stable, and the swing amplitude becomes small. The RMS temperature 
decreases at r = 0. This phenomenon was also reflected in the maximum RMS 
temperature, as shown in Fig. 2(f). When zf < 0.6, the maximum RMS temperature was 
observed at r*= 1. With increasing height, it moved to the center (r* = 0). Evidently, the 
RMS temperature is the highest at r*= 0 when zy = 0.6 in Fig. 2(e). When the position 


approaches the plate surface, the maximum RMS temperature was restored to 7’ = 1 at 
zf = 0.8. In the flow field, hot and cold fluid flowed along the positive axial direction. 
When the fluid reached the plate surface, the direction of fluid flow changed. The fluid 
flowed along the radial directions. Therefore, the junction of the cold and hot fluids 
moves outward. The maximum RMS temperature was observed at r* = 2 on the plate 
surface. Furthermore, the maximum RMS temperature at the position 0.2 < zf < 0.8 is 
much larger than that on the plate surface; this implies that the temperature-fluctuation 
attenuation is rapid near the plate surface. 


5. Effect of plate height on temperature fluctuation 


To experimentally investigate the effect of plate height, the plate was placed above 
the jet nozzle at different heights, and the temperature fluctuations on the lower surface 
were measured. The height range was 30-110 mm, denoted as Z3—Z9 in Fig. 1(e). 


5.1 Transient temperature changes with height 

Compared with zp = 50 mm, the cold fluid expanded and contracted with time at 
different plate heights. However, changes in the swing and shape of the jet were 
insignificant when the plate height zp = 30 mm. As the height increases, the jet swing 
and temperature distribution changes become increasingly pronounced, as shown in Fig. 
3. 

The transient temperature at zp = 30 mm and its frequency domain representation 
in the radial direction are shown in Fig. 4. The temperature fluctuation was the most 
severe at r* = 2 between the hot and cold fluids at this height. The fluctuation range was 
approximately 0.2—0.9. The frequency domain (Fig. 4(a)) representation shows that the 
temperature fluctuation has no prominent peak frequency. However, the temperature 
fluctuations were more severe at lower frequencies. At r= 2, the frequency of the 
temperature fluctuation was mainly distributed within 10 Hz. At r*= 0, the amplitude 
above 5 Hz is almost zero. 
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Fig. 3 Dimensionless transient temperature distribution on the plate surface 
The transient temperature variations at plate heights zp = 40, 50, 70, and 90 mm 
are shown in Fig. 4. The ranges of temperature fluctuations at different positions 
exhibited different changes in the axial direction. At the center (7 = 0), the temperature- 


fluctuation range changed from 0—0.1 at z»= 30 mm to 0. 0.7 at zp= 50 mm and 0.7— 
0.8 at z» = 90 mm. The temperature-fluctuation range at r’= 2 decreased from 0.2-0.9 
at zp = 30 mm to 0.7-0.8 at zp = 90 mm. The temperature fluctuation at r*= 5 was from 
0.9-1 at zp» = 30 mm to 0.75—0.85 at zp = 50 mm and finally reached 0.7—0.8 at zp = 90 
mm. During the entire process, the fluctuation range remained almost unchanged. Only 
the average temperature decreased slightly. The variation in the temperature-fluctuation 
range was the result of cold and hot fluids mixing and jet instability. The variation in 
temperature fluctuation with height was analyzed using the average temperature and 
temperature fluctuation intensity. 
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Fig. 4 Transient dimensionless temperature in time and frequency domains 


5.2 Average temperature changes with height 

The average temperature of the center gradually increased in the positive axial 
direction. The outside temperature gradually decreased, as shown in Fig. 5(a). The 
temperature distribution was almost circular. Therefore, the temperature in the radial 
direction was used for further analysis. 
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(b) intensity of temperature fluctuations 


Fig. 5 Dimensionless average temperature on the plate surface 


The dimensionless average temperature of the plate surface in the radial direction 
is shown in Fig. 6(a). The temperature in the radial direction was symmetrical. However, 
because of the non-coaxial nature of the testing bracket, its symmetry region was 
between r“ = 0 and r= -1. The average temperature at the center was low and increased 
gradually with the radius. It reached a constant value at r= 5. As the height increased, 
the temperature at the center increased. The temperature at each radial point was 
identical. This indicated that the mixture became uniform. 
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Fig. 6 Dimensionless average temperature and intensity of temperature fluctuations on the 
plate surface 

The dimensionless average temperature variation along the axial direction was as 
shown in Fig. 6(b). The temperature at 7" = 0 increased continuously with increasing 
height. The position at r* = 2 showed a relatively slow increase. However, the 
temperature at r= 6 decreases in the range of 30 mm < zp < 70 mm. When z, > 70, it 
remained almost constant. This demonstrates the effect of height on fluid mixing. At zp 
= 30 mm, the average temperature at 7’ = 0 is approximately 0. This implies that the 
cold fluid at the center did not start to mix. The temperature at each point changed 
rapidly within the range 30 mm < zp < 70 mm. This implies that the fluid mixes rapidly 
at this height. When z, > 70 mm, the temperature gradually varied. The temperatures at 


each point tended to be the same, indicating that the mixing process was complete. 


5.3 Temperature-fluctuation intensity changes with height 

The position of the maximum temperature-fluctuation intensity on the plate 
surface evidently changed with height. When z, = 30 mm, an annular shape was 
observed with a strong temperature fluctuation between the hot and cold fluids. The 
RMS temperature at the center was relatively low, as shown in Fig. 5(b). However, 
when zp = 40 mm, the RMS temperature above the nozzle showed little difference and 
only decreased significantly when r“> 2. As the height increased, the RMS temperature 
gradually decreased, and the temperature-fluctuation intensity at each point gradually 
became uniform. 

The RMS temperature variation in the radial direction are shown in Fig. 6(c). The 
points with higher temperature-fluctuation intensity at different heights are mainly 
distributed in the circular area above the nozzle in the range r’ < 4. The temperature- 
fluctuation intensity is determined by the jet instability and hot and cold fluid mixing 
with height. As the height increases, the amplitude of the jet sway increases. For low 
values of height, the point of the most severe temperature fluctuation is between the hot 
and cold fluids (r°= 2). However, when the height is higher, the intensity of temperature 
fluctuation in the range of r*< 4 tends to be consistent. However, with an increase in 
the axial direction, the jet mixing becomes more uniform, and the amplitude of the 
temperature fluctuation decreases gradually. Under the combined action of these two 
factors, the temperature-fluctuation intensity first increases and then decreases at r= 0, 
and it decreases continuously at 7 = 2. It remains almost unchanged when r’ = 6, as 
shown in Fig. 6(d). This implies that when the plate is located at different heights above 
the jet, the points affected by the temperature fluctuations that need to be analyzed are 
entirely different. When the height is low, the analysis must focus on the locations 
between the hot and cold-fluid nozzles. When the height is high, a series of positions 
from the center of the cold-fluid nozzle to the position of the hot fluid nozzle must be 
comprehensively considered. The analysis of the temperature-fluctuation range in 
Subsection 5.1 also highlights this phenomenon. 


6. Conclusions 


In this study, fluid and plate surface temperature fluctuations were experimentally 
investigated using a setup based on a coaxial-jet model. Essentially, the behavior of the 
temperature fluctuation of the fluid and the plate surface was studied. Furthermore, the 
variation in the axial direction was analyzed. The following conclusions were drawn 
from this study: 

(1) An experimental setup was constructed. The combined standard uncertainty of 
the temperature measurement is u-(7T) = 0.311 K, and the expanded uncertainties of the 
flow measurement Ucreis(V) = 0.42%-8.4%. The setup exhibited sufficient stability and 
accuracy to study the temperature fluctuations of the fluid and plate surface. 

(2) The temperature fluctuations were mainly caused by the instability of the jet. 
The instability of the fluid was primarily manifested in the expansion and contraction 


of the center of the cold fluid over time, the swing of the center position of the jet, and 
the shape change of the jet. These three factors caused temperature fluctuations on the 
plate surface. With an increase in the plate height, the jet swing and shape change 
became more pronounced. 

(3) The dimensionless average temperatures of the fluid and plate surfaces were 
symmetrical around the center of the cold fluid. For the fluid, it first increased and then 
decreased along the radial direction; further, it varied linearly along the axial direction 
within zy = 0.2-0.8. On the plate surface, it increased continuously as the radius 
increased. It gradually increased at the center with increasing plate height. The average 
temperature reflects the fluid mixing. Fluid mixes rapidly during 30 mm < z»< 70 mm 
and slowly when zp > 70 mm. 

(4) The maximum temperature-fluctuation intensity on the plate surface was much 
smaller than that of the fluid. The temperature-fluctuation attenuation was rapid near 
the surface of the plate. The intensity of the temperature fluctuation on the plate surface 
was affected by its radius and height. It first increased and then decreased at r* = 0, 
whereas it continuously decreased at r*= 2 and remained unchanged at r*= 6 as the 
plate height increased. When the plate height was low (zp < 30 mm), the analysis should 
focus on the locations between the hot- and cold fluid nozzles. When the height is high 
(Zp > 40 mm), a series of positions from the center of the cold-fluid nozzle to the hot 
fluid nozzle must be comprehensively considered. 
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Nomenclature 
CHT Conjugate heat transfer 
DSM Dynamic Smagorinsky model 
IHCP Inverse heat conduction problem 
LES Large eddy simulation 
RANS Reynolds-averaged Navier-Stokes 
RMS Root mean square 
SGS Subgrid-scale 
URANS Unsteady Reynolds-averaged Navier-Stokes 
AT Hox) Half-width of the data acquisition card, K 
AV powmeter) Half-width of electromagnetic flowmeter, m?/h 
r Radius, mm 
Cold Radius of the cold-fluid inlet, mm 
r“ Dimensionless radius r“= r 
T Temperature, K 


T cold Inlet temperature of the cold fluid, K 
THot Inlet temperature of the hot fluid, K 
T? Dimensionless temperature T”= ae 
THor-T Cold 
Dimensionless average temperature 
Ti IX 
oe T hig 3 A r 
n 
i=l 
Dimensionless root-mean-square temperature 
Th ix 2 
Tius = [=> (TT) 
n 
i=l 
t Time, s 
OT hermocouple) Expanded uncertainty of thermocouple, K 
u(THroKi) Standard uncertainty of data acquisition card, K 
u(T ihermocouple) Standard uncertainty of thermocouple, K 
u(V) Standard uncertainty of flow, m/h 
UV powmeter) Standard uncertainty of flowmeter, m3/h 
Combined standard uncertainty of temperature measurement, K 
u(T) | 2 2 
u-(T) = U(Tihermocouple) +u(T Hioki) 
Combined standard uncertainty of inlet temperature measurement, K 
; 2 2 
UclTin) = [(Tinermocouple) | +[u(T ricki) 
01) = eo =r: 
Relative expanded uncertainty range of flow 
Ucrels( V) us( y) Us( V) 
u =k ~ 
cres V) ( V nax V nin ) 
V flow, strip/h 
Zf Measurement height in the fluid, mm 
Zp Plate height, mm 
zy Dimensionless measurement height in the fluid z=% 
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